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ABSTRACT 

We report the detection of sub-Saturn-mass planet MOA-2008-BLG-310Lb 
and argue that it is the strongest candidate yet for a bulge planet. Deviations 
from the single-lens fit are smoothed out by finite-source effects and so are not 
immediately apparent from the light curve. Nevertheless, we find that a model 
in which the primary has a planetary companion is favored over the single-lens 
model by Ax^ ~ 880 for an additional three degrees of freedom. Detailed analysis 
yields a planet/star mass ratio q = (3.3 ± 0.3) x 10""^ and an angular separation 
between the planet and star within 10% of the angular Einstein radius. The 
small angular Einstein radius, 6e = 0.155 ± 0.011 mas, constrains the distance 
to the lens to be > 6.0 kpc if it is a star {M^ > 0.08 M©). This is the only 
microlensing exoplanet host discovered so far that must be in the bulge if it is a 
star. By analyzing VLT NACO adaptive optics images taken near the baseline of 
the event, we detect additional blended light that is aligned to within 130 mas of 
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the lensed source. This light is plausibly from the lens, but could also be due to 
a companion to lens or source, or possibly an unassociated star. If the blended 
light is indeed due to the lens, we can estimate the mass of the lens, = 0.67 ± 
0.14 Mq, planet mass m = 74 ± 17M^, and projected separation between the 
planet and host, 1.25 ± 0.10 AU, putting it right on the "snow line". If not, then 
the planet has lower mass, is closer to its host and is colder. To distinguish among 
these possibilities on reasonable timescales would require obtaining Hubble Space 
Telescope images almost immediately, before the source-lens relative motion of 
/i = 5masyr~^ causes them to separate substantially. 



1. Introduction 



Over the past 5 years, gravitational microlensing has led to the discovery of several 
exoplanets that would not be detectable by any other method currently available. Because 
it does not rely on light coming from the planet or the host star, microlensing is able to 
detect planets at several kpc, probing even into the center of the Galaxy. Thus microlensing 
has the potential to determine the demographics of planets orbiting hosts from two distinct 
stellar populations, bulge stars and disk stars, which is critical for understanding the Galactic 
distribution of planets and may allow one to constrain the time-history of planet formation 
in the universe. 

Standard models of the spatial and velocity distributions of stars in the Galaxy predict 
that roughly 2/3 of all microlensing events of st ars in the bulge arise from bulge lenses, 
i.e., stars in the bulge (IKiraga fc Paczyhskil 119941 ). In light of this, one might expect that 
planetary detections via microlensing would be more frequent in the bulge than in the 



disk. On the contrary, of the eight microlensing planets discovered so far (IBond et al. 
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Dong et al.ll2009bl ). five have measured or constrained lens distances, and four of these 



are clearly in the foreground disk while none are unambiguously in the bulge. The distance 
is not well constrained for the remaining three planets (OGLE-2005-BLG-169Lb, OGLE- 
2005-BLG-390Lb, and OGLE-2007-BLG-400Lb), and none have been definitively identified 
as bulge planets. The low detection rate of bulge planets may arise from a selection bias 
that favors the longer events that preferentially arise from disk lenses, or it may reflect the 
underlying Galactic distribution of planets. Here we present the analysis of a planetary sig- 
nature in microlensing event MOA-2008-BLG-310, the strongest candidate for a bulge planet 
to date. 



In Section [2l we discuss the observations and data reduction. In Section [3l we fit these 
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data to single-lens and planetary models. In Section H] we discuss our treatment of limb 
darkening of the source, which is important because it potentially affects the light curve 
at the times of maximum deviation due to the planet. We measure the Einstein radius in 
Section O and thereby constrain a combination of the lens mass and distance. To obtain 
a second, independent constraint on these two quanities, we first search, in Section [6l for 
"microlensing parallax" effects, but these prove too small to be observed. The analysis of 
images from SMARTS CTIO, VLT NACO, and IRSF, detailed in Section El does however 
reveal excess light aligned with the event. In Section [H] we discuss the four possible sources 
of this excess light: the lens, a companion to the lens, a companion to the source, or an 
ambient star. Finally, Section [9] describes how future observations with HST or adaptive 
optics (AO) may help characterize the host and its planet. 



2. Observations 

The Microlensing Observations in Astrophysics (MOA) collaboration detected microlens- 
ing event MOA-2008-BLG-310 [(RA,Dec) = (17:54:14.53,-34:46:40.99), (l,b) = (355.92,-4.56)] 
on 6 July 2008 (HJD' = HJD - 2450000 = 4654.458). MOA issued a high-magnification 
alert two days later, about 12 hours before the event peaked. The color and magnitude 
of the source indicate that it is a G type star i n the Galactic bulge, a result confirmed by 



high-resolution spectroscopy (ICohen et al.ll2009l ). 



The Microlensing Follow Up Network (/iFUN) began to intensively monitor this event 
at HJD' = 4656.026, less than 9 hours before the peak. The minimum predicted peak 
magnification was A^ax > 80 but the best-fit model at the time was consistent with formally 
infinite magnification, so the event was given high priority. Observations were taken by 
six observatories, MOA (New Zealand) 1.8 m /, /xFUN Auckland (New Zealand) 0.41m 
i?, /iFUN Bronberg (South Africa) 0.36 m unfihered, /iFUN SMARTS CTIO (Chile) 1.3 m 
/, V, H, MiNDSTEp La Silla (Chile) 1.54 m /, and PLANET Canopus (Tasmania) 1.0 m 
/. Only one observatory, /iFUN Bronberg, was positioned to see the peak of the event. 
Nevertheless, this observatory provided very complete coverage. /iFUN Bronberg took a 
total of 973 observations over the period 4656.21 < HJD' < 4656.55, recording the peak 
and all interesting anomalies. The high density of these observations allows us to bin the 
/iFUN Bronberg data without compromising the time resolution. The binned data points 
(as seen in Fig. [H below) occur every 2.5 minutes over the peak, whereas each planetary 
feature spans roughly an hour. /iFUN SMARTS took a total of 49 images in the /-band, 
275 images in the i/-band, and 6 in the V-band. The /iFUN SMARTS observations overlap 
/iFUN Bronberg by about 2 hours, starting after the peak and providing additional coverage 
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of the last planetary deviation. 



For /iFUN observatories, we primarily use difference imaging analysis (DIA) (jWozniak 



20001 ). but also use DoPhot reductions of /xFUN SMARTS if-band data to investigate the 
light that is blended with the sour ce. MOA data we re reduced using the standard MOA dif- 
ference imaging analysis pipeline (IBond et al.ll200ll ). PLANET Canopus data were reduced 
using the pySIS2 pipeline, based on the ISIS 2 code of lAlardI (120001 ). The error bars for all 
data are renormalized so that P^i' degree of freedom for the best-fit planetary model is 
close to unity. 

Being unfiltered, the /iFUN Bronberg data are subject to a differential extinction cor- 
rection because the source has a different color than the mean color of the reference frame 
used by DIA. We measure this effect from t he light curves of st able stars having the same 
color as the lens, and thereby remove it. See iDong et al.l (l2009al ). 



3. Microlens Model 

MOA-2008-BLG-310 was initially modeled as a single lens event. The single-lens model 
light curve fits the data reasonably well, showing pronounced finite-source effects in the 
rounding of the peak but no obvious anomalies. The event reached a maximum magnification 
^max ~ 400, making it a good candidate for planet detection although the finite-source 
effects work to smooth out any planetary deviations. Figure [1] shows the light curve and 
the residuals to the the best-fit single-lens and planetary models. The model allows us to 
roughly determine several parameters pertaining to the general structure of the light curve; 
^0; ""o; ^E, and p. Here, to is the time of minimum separation between the source and lens, Uq 
is the minimum separation in units of the Einstein radius, p is the radius of the source in 
the same units, and is the Einstein crossing time. We find that the source crossing time, 
= ptE, is better constrained than p, and so we report this parameter as well. 

A close look at the residuals from the single-lens fit (middle panel of Fig.[T]) reveals signif- 
icant structure indicating that the underlying lens model is more complicated than a simple 
single lens. In particular, short timescale deviations near the peak of high-magnification 
events are typically caused by a planetary or binary companion. In these cases, the caustic 
structure is extended, as opposed to the simple point in the case of an isolated lens, leading 
to deviations from the single-lens form as the source crosses the caustic. The most impor- 
tant features, a short spike in the residuals just before the peak (HJD' 4656.34), and a short 
dip just after (HJD' 4656.48) are completely covered by unfiltered observations from /iFUN 
Bronberg. The second of these features is confirmed in /-band data from /xFUN SMARTS. 
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SMARTS i7-band observations qualitatively show the same deviation despite suffering from 
larger scatter. Since the higher quality /-band data cover the same portion of the light 
curve, i/-band data are not used in the derivation of model parameters. La Silla J-band 
data further confirm the last half of the second feature. Because Bronberg provides the most 
crucial coverage of the anom alies, we conduct thre e additional independent reductions of 
Br onberg data using DoPhot (ISchechter et al.lll993l ). the DIA reduction pack age de v elope d 
bv bond et aP J200lh . and the phSIS2 pipeline, based on the ISIS 2 code of lAlardI J200oh . 
All three confirm the structure of the light curve in this critical region. 

Note that the pronounced misalignment of the Bronberg and yuFUN SMARTS data in 
the middle panel is real: because /xFUN SMARTS data do not cover the peak, the fs and 
/b parameters are permitted much more freedom to match the single-lens model than is the 
case for Bronberg. 

The relatively low amplitude of the residuals from a single-lens model, along with the 
fact that these residuals are apparent over most of the duration of the source diameter 
crossing, generally indicate that the central caustic str ucture due to a companion to the len s 
is only magnifying a fraction of the source at one time (IGriest fc Safizadehlll998l : lHanll2007l ). 
This suggests that w, the "short diameter" or "w idth" of the centra l caustic is smaller than 
or comparable to the diameter of the source (see lChung et al.ll2005l ). Prominent deviations 
from the single-lens model occur where the limb of the source enters and exits the caustic. 



This b ehavior, which is qualitatively very similar to that of MOA-2007-BLG-400 (jPong et al. 



2009bl ). prompts us to investigate possible two-point-mass lens (planetary or binary) models. 



3.1. Searching a Grid of Lens Geometries 



The finite-source two-point-mass lens magnification calculations are carried out using the 
improved magnification map technique of iDong et al.l (l2006l . l2009bl ). w hich is optimized fo r 
high-magnification events. The fitting procedure follows closely that of iDong et al.l (l2009bl ). 
The initial search for two-point- mass lens solutions is conducted over a grid of three pa- 
rameters: the short-caustic width w, companion/star mass ratio q, and the angle of the 
source trajectory relative to the companion/star axis a. Since w is a function of q and 
the companion/star separation d, this is equivalent to also fixing d at various values. The 
remaining parameters {tQ,UQ,t^,p) are allowed to vary. Two additional parameters for limb 
darkening are given fixed values, as will be discussed in Section HI The source flux fs and 
the blended flux fb are fit independently for each filter and telescope. We use Monte Carlo 
Markov Chain (MCMC) to minimize with respect to {to,uo,t-E,p) at each of the {w, q, a) 
grid points. There is a well known degeneracy such that for g -C 1, planet /star separations 
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d and will produce almost identical central caustic structures and consequently indistin - 
guishable light curves for high magnification events such as this (iGriest &: Safizadehlll998l ). 
We explore a {w,q) grid for each geometry, searching the d > 1 (in units of the Einstein 
radius) regime for 'wide' solutions and the d <1 regime for 'close' solutions. 



3.2. Best-Fit Model 

An initial search for two-point- mass lens solutions is conducted over the range of caustic 
widths —3.5 < \ogw < —1.0 (in units of the Einstein radius), companion mass ratios —5.0 < 
log 5 < 0, and source trajectory angles < a < 27r in the two separate regimes d>l and 
d < 1. This initial search gives us fairly good estimates of the best-fit parameters and the 
location of the minima in terms of w and q (and hence also d). For this particular event, 
however, w and q turn out to be highly correlated. We conduct a refined search over a grid 
in (c?, g) instead of {w,q), and we also allow a to vary as a MCMC variable, rather than 
discretely. The solid black lines in Figure [2] show Ax^ = 1,4,9 contours in the {d,q) plane 
for the wide (top) and close {bottom) solutions, respectively. For the wide solution, the 
minimum occurs at ci = 1.085 ± 0.003 and q = (3.31 ± 0.26) x 10"'^. The close solution 
minimum occurs ai d = 0.927 ±0.003 and q = (3.20 ±0.26) x 10^^. The mass ratio indicates 
that the companion to the lens is in fact a planet. As expected, we recover the d ^ d~^ 
degeneracy. The wide solution is favored by just Ax^ = 2.06, indicating that the wide/close 
degeneracy cannot be clearly resolved in this case. The best-fit parameters for both wide and 
close solutions are recorded in Table [H Two independent algorithms were used to explore 
parameter space and both returned essentially identical best fits. 

The wide and close planetary models qualitatively explain several features of the single- 
lens residuals. The lower panel of Figure [3] shows the extended source at key points in 

Table 1. MOA-2008-BLG-310 Best-Fit Planetary Model Parameters 



d g X 10^ to(HJD') uq x 10^ tE(days) a(rad) t*(days) X^/dot 

Wide 1.085 3.31 4656.39975 3.00 11.14 1.21 0.05487 2891.40/3050 

Close 0.927 3.20 4656.39975 3.01 11.08 1.21 0.05483 2893.46/3050 

Error 0.003 0.26 0.00005 0.14 0.50 0.02 0.00009 - 



single-lens x^/dof = 3773.23/3053 
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time on its trajectory. The nearly identical central caustics generated by the best-fit wide 
and close models are both shown. The most prominent features in the residuals, shown 
in the upper panel of Figure El occur as the limb of the source crosses the caustic. The 
positive and negative spikes that are most evident from the raw data (features 2 and 4 
of Fig. [3]) coincide with the limb of the source entering the strong curved portion of the 
caustic and exiting the weaker straight segment. Residual patterns like this, characterized 
by short duration perturbations separated by a relativ ely flat region, are typical of planetary 



lens systems affected by strong finite source effects (iGriest fc Safizadehl 1 19981 : iDong et al. 



2009bl : iHan &: Kinj|2009l ). The bottom panel of Figure [T] shows the residuals to the best- 



fit wide planetary model. The deviations from the point-lens model that initially indicated 
that the lens was not being accurately modeled are no longer apparent. The planetary model 
decreases by ~ 880 for three additional degrees of freedom, making this a strong detection 
of a Saturn mass-ratio planet/star system. 



4. Limb Darkening 



Since the primary deviations in the light curve occur at the limb of the star, it is 
important to examine the effects of limb darkening on the planetary solution. For this 
purpose we adopt a surface brightness profile of the form 



So 



1 — -{cos'd) 



A 



l_^(cosl/2^) 



where 'd is the angle between the normal to the surface of the star and the line of sight 
( lAn et al.ll2002l ). This is somewhat more complicated than the typically used linear limb 
darkening, but it is justified by an improvement in the fit of Ax^ ~ 10. We simply fix the 
MOA, CTIO, and La Silla limb-da rkening parameters at (F, A) = (0.077,0.549) correspond- 
ing to {c,d) = (0.099,0.584) from IClaretl (120001 ). These parameters pertain to a star with 
Teff = 5750 K and logg = 4.0, i.e., a post-turnoff G star, corresponding to the (y — I)o = 0.69 
and Mj = 3.46 that we derive from the color-magnitude diagram by assuming that the source 
suffers the same extinction and is at the same distance as the bulge clump. 

However, because /iFUN Bronberg provides the bulk of the observations covering the 
peak of the event and the deviations at the limb of the source, it is most critical that the limb 
darkening be accurately modeled for these data. We first determine the effective bandpass 
of Bronberg (which is unfiltered) by making a color-color diagram of stars in the field with 
colors similar to that of the microlensed source. We find that 



A(i?Bron-/) = 0.50A(l^-/), 



(2) 



- 10 - 



i.e., almost exactly what would be expected for standard R band. We ther efore begin by 
adopting (r,A) = (0.166,0.543) corresponding to {c,d) = (0.204,0.557) from Iciaretl J200oh 
for R band. As a check on this procedure, we also allow F and A for Bronberg data to vary 
(along with most other pa rameters), but still holding the /-band limb-darkening parameters 
fixed at the I Clare tl (120001 ) values for the other observatories. The best-fit models for wide 
and close planet/star separations have (F, A) = (-0.200,1.277) and (r,A) = (0.069,0.732) 
respec tively. We find that the resulting surface brightness profiles are similar to those defined 
by the I Clare tl (120001 ) parameters. 



We further investigate the effect of limb- darkening parameters on the final results by 
comparing likelihood contours for {d,q) for the two cases just described. The Ay^ contours 
for the close and wide planetary models with limb darkening fixed at the IClaretl (l2000l ) 
values are shown in Figure [2] as the solid lines. These contours are similar to the dotted 
lines in Figure [2] generated by allowing the parameters for limb darkening to vary freely. 
Most importantly, the best-fit values of d and q change by much less than one sigma. This 
justifies fixing the parameters for limb darkening at the I Claret! ( l2000l ) values for all models 
that follow. 

The temperatu re and metallicity we use are slightly different than those obtained by 
Cohen et al.l (120091 ) from spectroscopy of the event. The spectroscopy yields (F,A) = 
(0.105,0.535) ((c,d) = (0.133,0.564)) for / band and (F,A) = (0.203,0.517) ((c,rf) = 
(0.248,0.525)) for R band. Inserting these values into the wide and close models, we ob- 
tain best-fit parameters well within one sigma of those recorded in Table [T] and contours 
essentially identical to those in Figure [2l 



5. Measurement of Angular Einstein Radius 9e 



The color and magnitude of the source allow us to determine its angular radius which 
in turn can be used to place constraints on the lens mass and lens-source relative parallax. 
We begin by measuring the color and magnitude of the source and the clump centroid in the 
calibrated CTIO field, [{V - I), /],ource = (1-48 ± 0.01, 19.28 ± 0.05) and [{V - I), J]ciump = 
(1.84, 15.62). The source magnitude is derived from the microlens model and is the same for 
the close and wide solutions. At Galactic longitude / = —4.08, the angle of the bar-shaped 
bulge places the peak of the red clump density behind the Galactic center by 0.05 mag 
( iNishiyama et al.ll2005l ). Assuming a distance to the Galactic center of 8 kpc, the dereddened 
position of the clump is then [{V — /)o, /o]ciump = (1-05, 14.37). Thus the extinction toward 
the source is [E{V — I),Ai] = (0.79,1.25). We find the dereddened col or and ma g nitud e 
of the source, [{V - /)o, /o] source = (0.69, 18.03). Applying the method of IYoo et all tooi ). 
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we convert {V — I) to {V — K) using the color-color relations of iBessell &: BrettI (119881 ) . 
and we obt ain \{V — K)n, i^njsnnr rp = (1-48, 17.24). We then use the color /surf ace-brightness 
relations of iKervella et al.l (120041 ) to calculate the angular source radius, 



= 0.76 ± 0.05 /ias, 



(3) 



which (as with the next three equations) applies equally to both the wide and close solutions. 
The source crossing time is 



U = ptE = 0.05485 ± 0.00009 days, 
which implies that the (geocentric) proper motion, /igeo = 6**/^*, is 

Pgeo = 5.1 ± 0.3 mas yr^"^ 
The inferred Einstein radius, 6'e = /^geo^E, is then 

9e = 0.155 ±0.011 mas 



(4) 



(5) 



(6) 



The fractional uncertainties in 9^, 6e, and comparable, a typical result for point-lens 



events with finite source effects (lYee et al. 
source-lens relative parallax vTrci (see 



Gould 



20091). We can relate the lens mass to the 



2000 for details). 



Mr 



(7) 



where k = AG/c^ AU ~ 8.1 mas Mq^ If we require that Ml > 0.08 M© (that is, if the lens 
is a star) then it follows that vTrei < 37;uas. Assuming Ds > 8kpc (as discussed above, the 
Galactic bar at / ~ —4° lies behind the Galactic center), this gives a lower limit on the 
distance to the lens > 6kpc. We conclude that if the lens mass is above the hydrogen 
burning limit, then it must be located in the Galactic bulge. In order to verify the bulge 
location of the lens, we would need another independent relation between the lens mass and 
distance. This could be obtained by measuring either the microlensing parallax or the flux 
from the lens. 



6. Parallax 



Determining the microlensing parallax tte gives us a ind ependent relationship between 
the lens mass and source- lens relative parallax (lGouldll2000l ). The magnitude of the vector 
is given by. 
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while the direction is the same as that of ^tgeo; the lens-source relative proper motion in 
the geocentric frame. In combination with the independent relation between Ml and vTrei 
obtained from the proper motion of the source, it would be possible to give physical values to 
both of these parameters. With this goal in mind we examine the effects on the light curve 
from two sources of parallax. Orbital parallax is caused by the acceleration of the Earth 
on its orbit. Terrestrial parallax arises from two or more widely separated observatories 
simultaneously observing a slightly different light curve due to their different vantage points. 
For this event orbital parallax is not expected to be detectable since the timescale is so short 
(tE = 11.1 days). We expect terrestrial parallax to be poorly constrained as well. Earth- 
based parallax measurements require t hat short duration cau stic crossings be observed by 
two or more telescopes simultaneously (iHardy &: Walkerlll995l ). While /xFUN Bronberg and 
/iFUN SMARTS both observed the second prominent deviation as the limb of the source 
exited the caustic, this feature is washed out by finite source effects. We once again search the 
{d, q) grid, allowing the north and east components of both terrestrial and orbital parallax 
to vary as additional MCMC parameters. We also test the case of the source-le ns minimum 



separ ation uq —Uq as this is a known degeneracy in determining parallax (jSmith et al. 



20031 ) . For the four cases (±mo, close/wide), the reduction in ranges from 2 to 6, i.e., 
barely different from the Ax^ = 2 expected from reducing the degrees of freedom by 2. 
The marginal detection of parallax at Ax^ = 6 favors tte — 4. Such a large parallax yields 
TTrei = tte^e = 0.65 mas and lens mass Ml = 0.005 M0. We do not give much weight to 
this marginal parallax detection and the free-floating planet solution (with sub-Earth mass 
moon) it implies since from previous experience we ha ve found that such sma ll Ax^ could 
easily be produced by low- level systematics (see also iPoindexter et al.l 120051 ) . Hence we 
obtain essentially no new information, and, as our results are consistent with zero orbital 
and terrestrial parallax, we set vte = except where explicitly indicated. 



7. Blended Light 



From the best-fit model we obtain a measure of how much light is being lensed in the 
event, in other words the flux of the unmagnified source. In addition to the source flux, there 
is blended light that is not being lensed, which may come from unrelated stars along the line 
of sight, companions to the source or lens, or the lens itself. An alternate route t o obtaining 



the lens mass and d istance is possible if the flux from the lens can be isolated (IHanl 12005 
Bennett et aPboOTh . 



We have if-band images of the event taken from CTIO. We have additional post-event 
JHK Infrared images of MOA-2008-BLG-310 taken with the IRSF telescope in South Africa 
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on 2008 August 4 and the adaptive optics system NACO on the ESO VLT on 2008 July 
28 under ESO Program ID 081.C-0429(A). The pixel scales are respectively 0.27, 0.45 and 
0.027 arcsec. A log of the VLT NACO and IRSF observations is given in Table[21 As detailed 
in Appendix [Bl IRSF serves as a bridge between NACO and CTIO, being wider than the 
former and deeper than the latter. 

The NACO image reveals two additional stars in the vicinity of the source that are 
unresolved by the observations used in the light curve analysis. One of these is 3 mag 
brighter than the slightly magnified [A = 1.09) source and 0.85 arcsec away (star 3 in Fig. [S]) 
while the other is 0.2 mag brighter and 0.5 arcsec away. To definitively identify the source 
from among this group, we create a template image from the best CTIO /-band images and 
subtract this template from 20 different astrometrically aligned, good-seeing images near 
the peak of the event. The magnified light of the source is isolated on the subtracted image 
because the contribution from other stars is removed. Thus the relative astrometry of the 
source is very precisely determined. DoPhot is used to find the positions of other stars on the 
template CTIO / and median NACO H images. We select 14 isolated stars common to both 
images and calculate the coordinate transformation from CTIO to NACO. The position of 
the source transformed to NACO coordinates is 13 ± 25 mas from the centroid of the target 
in Figure [51 The nearest neighboring star is 400 mas from the source position, and thus the 
identification of the source with the target on the NACO frame is very secure. 

We reduce the IRSF images following standard procedures, and measure the fluxes and 
positions of stars using the DoPhot software. The reduction of the NACO images is a more 
complicated procedure, and is detailed in Appendix |Al Our goal is to put NACO photometry 
of the target (blend + magnified source) on the CTIO photometric system, so that it can 
be compared with the source-only iJ-band flux, which is well-measured from the CTIO H- 



Table 2. Log of Observations 



Image 


Date 


hour 


FWHM 


JlRSF 


200J 


^-08-04 


18:05:25 


1.4 arcsec 


HlRSF 


200J 


^-08-04 


18:05:25 


1.4 arcsec 


-f^IRSF 


200J 


^-08-04 


18:05:25 


1.3 arcsec 


-^NACO 


200J 


^-07-28 


01:29:20 


0.15 arcsec 


-f^NACO 


200J 


^-07-28 


02:18:29 


0.13 arcsec 


-f^NACO 


200J 


^-07-28 


00:36:17 


0.15 arcsec 
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band light curve. In principle, this could be done using comparison stars common to NACO 
and CTIO H band. However, there are only two such stars, and they have relatively large 
photometric errors in CTIO photometry. Instead, we use a large number of common stars 
to photometrically align the CTIO and IRSF systems, which can therefore be done very 
accurately. We then align the NACO and IRSF systems based on 4 common stars, which 
have much smaller errors and consequently show smaller scatter than the CTIO stars. We 
align the IRSF system with 2MASS, allowing us to determine the calibrated magnitudes of 
the target and reference stars recorded in Table O The photometric calibration of IRSF, 
NACO, and CTIO images is discussed in greater detail in Appendix [Bl We stress that in the 
following discussion (as well as the Appendices), we deliberately work in uncalibrated -ffcTio 
magnitudes, since measurement of the fraction of blended flux depends only on relative 
photometry. We apply the calibration only at the end of this procedure to avoid introducing 
additional uncertainty into the fairly subtle differential measurements. 



7.1. Estimation of the target flux in H CTIO 

In the following, the term "source" refers to the star that was microlensed, while the 
term "target" refers to all the light that is aligned with the source in the NACO images. 
We calibrate the NACO if-band magnitude of the target via the route IRSF-to-2MASS and 
flnd i^NAco,caiib = 17.47 ± 0.05. Then using the IRSF-to-CTIO transformation, we convert 
the measured NACO flux into the instrumental CTIO system, -fftarget,CTio = 21.29 ± 0.05. 
We stress that this indirect road NACO-to-IRSF-to-CTIO is actually the most accurate one 
to estimate the magnitude in the instrumental CTIO system. 

We also carry out the following independent check. We measure aperture fluxes fi on 
the NACO image for stars 1, 2, 3, and the target listed in Table |3l For the target, we 
correct the result for contaminating flux from star 3 and from another much fainter nearby 



Table 3. Photometric data for H CTIO, JHK IRSF, JHK NACO 



Star ID 


-f^CTIO 


-f^IRSF.calib 


-f^NACO.calib 


</lRSF,calib 


<^NACO,calib 


-f^IRSF.calib 


-f^NACO.calib 


1 


16.95 


13.094 


13.106 


13.872 


13.855 


12.88 


12.898 


2 


17.69 


13.834 


13.826 


14.230 


14.225 


13.77 


13.76 


3 




14.340 


14.352 


14.83 


14.884 


14.24 


14.246 


target 






17.47 




18.068 




17.349 
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star. Then using the same stars i = 1, 2, 3 on the IRSF image, we obtain an estimate 
of the target magnitude on the IRSF system: -fftarget;i,iRSF = -f^i,iRSF + 2.5 log (/i/ /target)- 
We take the average of these three estimates (whose standard error of the mean is only 
0.012 mag), and then apply the previously derived conversion from IRSF to CTIO. We find 

-f^target,CTIO = 21.27. 

As a further sanity check, we apply a similar procedure to compare the NACO and 
CTIO images directly. As stated at the outset, we expect that this will be less accurate both 
because there are only two viable comparison stars (1 and 2) and because the CTIO flux 
measurements are less accurate than those of IRSF. Nevertheless, we find a similar result: 
-f^target,CTio = 21.32, although with substantially worse precision. 

We finally adopt -fftarget.cTio = 21.28 ± 0.05, where the error bar reflects our estimate 
of the systematic error. Clearly, our two primary methods of estimating this quantity agree 
much more closely than this, but there still could be systematic effects common to both. We 
regard 0.05 mag as a conservative overestimate of the error. 

Inserting the if-band CTIO observations into the planetary model, we obtain the un- 
magnifled source flux, -ffsource,CTio = 21.55 ± 0.05 on the instrumental CTIO system. The 
error bar accounts for the the uncertainty in the flt by allowing all parameters, including 
parallax, to vary freely. However, for the purpose of determining the blend on the NACO 
image, we are more interested in the magnifled flux from the source at tNAco; the time the 
image was taken. The magnification is determined by the separation between the source and 
lens at tNAco, mnaco = (^naco — to)/tE in units of the Einstein radius. The unmagnified 
source fiux is anti-correlated with the Einstein crossing time ^e, so the dispersion in the mag- 
nified fiux is slightly smaller than the dispersion in the unmagnified fiux. The uncertainty 
in the magnified fiux is related to the model uncertainty in the unmagnified fiux (/g) by 

^ = 1 + ^ (9) 

which in the point-lens approximation (generally valid on the wings on the light curve) 
translates tc0 

^^^^«) 3-2A''[l-il-A~Y% (10) 



KM 



In our case, the analytical result '^^/f] = 0.77 is very close the result calculated using MCMC 
data. 



^In the limit of large A: 1 + ^^"^^ — -^^{^ + + -^^ji + •••)' which converges very quickly, even for 
A -2. 



- 16 - 



Figure [6] shows probability distributions for the magnified source fiux constructed from 
the MCMC chains. Without parallax (black histogram), the close and wide solutions give 
the same source magnitude at the time of the NACO image, -ffmagnified.cTio = 21.45 ± 0.04, 
while the best-fit solution with unconstrained parallax (gray) gives a fiux ~ 2% brighter. As 
noted in Section [6l the best-fit parallax imphes a planetary lens mass and is likely a spurious 
detection. If we constrain parallax so that the lens mass is at least 0.08 M0 (tte < 0.25), the 
best-fit magnified source fiux is identical to the case of no parallax. Thus our best estimate of 
the fiux strictly from the source is -ffmagmfied,CTio = 21.45 ±0.04, while the light aligned with 
the event on the NACO image is iftargct.cTio = 21.28±0.05. We consider 21.28 + 0.05 = 21.33 
to be a robust lower limit on the amount of light detected in the NACO image. Therefore 
excess light unrelated to the source is detected at the 3-sigma level. 



8. Constraints on the Origin of the Blended Light 

There are four possible causes of the excess fiux detected in the VLT NACO images: 
the lens, a companion to the lens, a companion to the source, or an ambient star unrelated 
to the event. We can estimate the probabilities of two of these possibilities (a companion 
to the source and an ambient star) fairly robustly using relatively well-constrained priors. 
However, estimates of the probabilities of the other two possibilites depend on relatively 
poorly-constrained assumptions about the distribution of planets orbiting hosts of various 
masses, including low-mass stars and brown dwarfs. Therefore, the relative probabilities of 
these four possibilities also depend on these assumptions, and so we cannot robustly distin- 
guish between them. We will therefore discuss the constraints we have on each possibility in 
turn, but will not attempt to assess which is the most likely interpretation. Rather, we will 
simply review the observations that are required to distinguish between these possibilities in 
Section 19. 3[ 

One constraint that applies to all four possibilities is that the excess light must lie within 
1 FWHM of the source. We find, by adding 17% of the target fiux at various positions, that 
separations greater than 1 FWHM would permit resolution of the excess fiux. The PSF of 
the VLT NACO image is slightly elongated: (140 mas x 124 mas), so we adopt 132 mas for 
present purposes. 
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8.1. Ambient Star 

From direct examination of the NACO H image, we find the density of stars within 
0.5 mag of the if-band magnitude of the excess hght to be 0.94 arcsec"^. Hence the prior 
probabihty that such a random star hes buried under the NACO image is 5.1%. 



8.2. Companion to the Source 



The source is a main-sequence G star. From Table 7 of iDuquennoy fc Mayorl (jl99ll ). 
we find that 9.4% of their sample (of 164 stars) have companions within the mass-ratio 
range of 0.57 < q < 0.76, i.e., the mas s range corresponding the ± 0.5 mag of the observed 
excess flux. However, from Figure 5 of IDuquennoy fc Mayorl (Il99ll ). 22% of companions lie 
outside ~ 1000 AU, the size of the NACO PSF projected on the source plane. A further 
3% have orbits shorter than ~ 3 days, which would have given rise to observable "xallarap" 
signals in the microlensing light curve. Hence, if bulge G stars are like the local sample, 
0.094 X 0.75 = 7.1% of them have companions within 0.5 mag of the observed excess flux, 
and lie at separations where they would not have been detected. This is comparable to the 
corresponding value for ambient stars. 



8.3. Companion to the Lens 

If the lens had a companion, it would induce shear on the lens's gravitational field, which 
would generate a small Chang-Refsdal (1979, 1984) caustic at the center of magnification of 
the lensing system. This would in turn produce spikes in the light curve at HJD' 4656.36 and 
4656.44, when the lens center-of-magnification crosses the limb of the source. The residuals 
to the planetary model (Fig. [1]) strongly limit any such spikes. To put this constraint on a 
quantitative basis, we fit the light curve to models that have two additional parameters, 0, 
the angle between the planet axis and the binary-companion axis, and if com, the width of 
the caustic induced by the companion, 

. Q'com / 1 1 \ 

Wcom = 4-^ , (11) 

^com 

where (gcom, c^com) are the mass ratio and separation of the companion. We hold (ifcom, 0) at 
a grid of fixed values and minimize with respect to all other parameters!^ 



^In practice, qcom is held fixed at 1, since we are probing the Chang-Refsdal (1979, 1984) limit, for which 
all caustics with the same w but various q are essentially the same. 
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This search reveals an improvement of Ax^ = — 7.3 for two additional degrees of freedom, 
with a best fit of (log ifcom; 0) = (—3.28,40°). This improvement is too small to claim a 
detection, since it could occur by chance with probability 2.6%, and could also be due to 
low-level systematics. Thus, while this test raises the tantalizing possibility that the excess 
light is due to a companion to the lens, we mainly focus on the 3 a upper limits to shear from 
a companion with Wcom < 1 x 10~^ over almost all angles. See Figure [71 By comparison, 

W^planet = 5 X 10"^ (sCC Fig. [3]). 



Equation (ITTI) can also be written 

4KMcom7rrel 



Wr 



where Mcom is the companion mass, and A^^com is its angular separation. If the lens is in the 
bulge and the excess light is due to its companion, then (see Section [H3D Mcom ~ 0.6 Mq, 
and so Wcom < 1 x 10^'^ implies 

A9c„m > 16 mas — = 16 mas f ^ . (12) 

37/ias VO.O8M0/ ^ ' 

For foreground- disk lenses, the limit on A^^com continues to grow, but much more slowly, to 
70 mas for (M, vTrd) = (0.003 Mq, 1 mas). 

Hence, in contrast to source companions, which are permitted over 4 decades of separa- 
tion, companions to the lens are restricted to 1-2 decades for bulge lenses and a somewhat 
narrower range for disk lenses. 



8.4. Lens Mass Estimate 

We stress that the excess light could be due to any of the three options above. The 
excess light could also be due to the lens. As we show below, this requires the lens star to be 
relatively massive and so quite close to the source. Such small source-lens distances are gen- 
erally disfavored by phase space and kinematical factors. However, as mentioned previously, 
evaluating the prior probability of this scenario requires adopting a specific assumption of 
frequency of planets as a function of host mass, which is poorly constrained. 

Under the assumption that the NACO blended light is due to the lens, we can estimate 
the lens mass using its inferred instrumental flux -ffiens.cTio = 23.38 ± 0.41, the measured 
Einstein radius 9e = -s/kMlvt^ = 0.155 ± 0.011, and an assumed range of mass-luminosity 
relations, which depend on the lens's unknown metallicity. While our actual calculation 
is fully self-consistent, the basic result can be understood intuitively as follows. For any 
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possible mass consistent with the observed flux, the lens-source distance will be quite small, 
^Ls ~ TTreil^g/AU = 300 pc (Ml/0. 7 Mq)"! relative to Dg, and the range of possible values 
due to different candidate masses is even smaller. Hence, we can just adopt Dls = 300 pc. 
For fixed Dls the lens lies, on average, -Dls/2 or 0.04 mag in front of the local bulge density 
peak (defined by the clump giants). At Galactic latitude b = —4.56, the dispersion of distance 
moduli of bulge sources is (5/ In 10) sin 6/0.6 = 0.29 mag for an adopted bulge flattening of 
0.6. We adopt an absolute magnitude for the clump of Mj:^ clump = —1.41 ±0.05, and so from 
the observed -ffciump,CTio = 17.35, we infer 

MH,lens = MH,clump + (^^lens.CTIO - i^clump,CTIo) + (0.04 ± 0.29) = 4.66 ± 0.50 



Then using six isochrones generated by the Dartmouth Stellar Evolution Database (IDotter et al 



20081 ) with [Fe/H] ranging from —0.5 to 0.5 and age ranging from 5Gyr to lOGyr, we esti- 



mate the lens mass to be Ml = 0.67 ± 0.14 M©. 

Assuming the excess light is indeed due to the lens and using the resulting estimate of the 
lens mass, we can now estimate the properties of the planetary companion. The planet mass 
is rrip = 74 ± 17M^, roughly 80% the mass of Saturn. Taking account of the uncertainties 
in both the distance to the lens and the angular Einstein radius, as well as the wide/close 
degeneracy, the projected separation between the planet and host star is 1.25 ± 0.10 AU. 

Note that if the blended light is not due to the lens, then the lens must be fainter than 
this light and so (unless it is a remnant) also of lower mass. The planet would then be of 
proportionately lower mass as well. 

9. Discussion 

9.1. Sub-Saturn Mass Planet — Candidate Bulge Planet 

Microlensing event MOA-2008-BLG-310 is one of only two published high magnification 
events to date for which the source is as large or lar ger than the centra l caustic. It bears 



many similarities to the other, MOA-2007-BLG-400 fiDong et al.ll2009bl ). Like that earlier 
event, the planetary perturbations in the light curve are not immediately apparent, having 
been smoothed out by finite-source effects. We find that a Saturn mass ratio planet/star 
model is nevertheless favored over the single-lens model by a significant reduction in 
(Ax^ ~ 880 for an additional three degrees of freedom). Using VLT NACO (together with 
IRSF) photometry, we definitively detect excess light blended with the source that is due to 
the lens, or a companion to the lens or the source, or an unassociated star. Regardless of 
the origin of this excess light, however, it places an upper limit on the lens flux and so on its 
mass. The planet's Saturn-like mass ratio therefore implies that it has a sub-Saturn mass. 



Our measurement of the angular Einstein radius 6e constrains a combination of the lens 
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mass and distance. We thereby conclude that if the lens is a star, then it must be in the 
bulge. 

We are not able to resolve the close/wide degeneracy in the geometry of the planetary 
system. However the separate solutions for the lens/star separation d differ only by a factor 
of 1.17. The d ^ 1/d degeneracy is not as severe in this case because the planet is located 
very close to the Einstein radius. 

Figure [S] shows the mass versus equilibrium temperature for planets that have been 
detected orbiting main-sequence stars via radial velocity, transits, direct imaging, astrometry, 
and microlensing. The position of MOA-2008-BLG-310Lb is shown under two assumptions. 
The hexagon symbol indicates its position assuming that the excess flux is due to the lens. 
We then obtain a host mass of Ml = 0.67±0.14Mo, and so a planet mass of 0.23±0.05 Mjup. 
The "tail" extending toward lower masses and colder temperatures assumes that the excess 
light is not from the lens (and so is due to a companion to the source or lens). The path of 
this tail is determined by the measurement of 6e = 0.155 mas, which constrains the product 
of the lens mass and lens-source relative parallax to be 6*1/ k = MTTrei = 3 Mq /las. 

9.2. Importance of Further Characterizing the Planet 

All possibilities for the lens identification are interesting. In particular, the host might 
be a brown dwarf or free-floating planet in the foreground disk. In this case, the Saturn- 
mass-ratio companion would actually be a moon. Future observations with the Hubble Space 
Telescope (HST) or adaptive optics (AO) could distinguish among these three possibilities. 
First, of course, mere detection of the light from the host would confirm it is a star and 
therefore that this is indeed a bulge planetary system, the first unambiguous such detection. 

If the excess light is due to the lens, then Figure [8] shows that this detection is probing 
a new part of parameter space. The previous eight planets detected via microlensing span 
a relatively wide range of mass from a few Earth masses to several Jupiter masses, but are 
largely located at relatively cold equilibrium temperatures of ~ 40—80 K, similar to the outer 
planets of our solar system. In contrast, the radial velocity and transit methods are generally 
only sensitive to sub-Saturn mass planets with relatively warm equilibrium temperatures of 
> 300 K. Therefore, little is currently known about the demographics of sub-Saturn mass 
planets with equilibrium temperatures between 100 — 300 K. MOA-2008-BLG-310Lb would 
make this the first microlensing planet to fall on the Snow Line. 

Of course, if the excess light is not due to the lens, then we have no hard information 
on the lens mass. However, the "tail" in Figure [H] indicates an interesting possibility: that 
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MOA-2008-BLG-310LB is a "cold Neptune" orbiting a low mass star, similar to several other 
such detections. 



9.3. Planet Characterization using HST or AO 

What are the prospects for characterizing the host and its planet? As we discussed 
in Section [6l the event contains essentially no parallax information. Hence, the only path 
toward measuring the lens mass and distance is direct detection of the host (or possibly its 
companion). For either the wide or close solution, the geocentric proper motion is /igeo = 
^e/^e = S.lmasyr"^. The heliocentric and geocentric proper motions differ by 

It^hei - Mgcol = |ve,±Frci = = 0.018 — masyr (13) 

where v^^± = 28kms~^ is the velocity of the Earth projected on the plane of the sky at the 
peak of the event. Hence, if the lens is luminous (M > 0.08 Mq), then the heliocentric and 
geocentric proper motions are essentially identical, and so the magnitude of the heliocentric 
lens-source relative proper motion is well determined. 

This known proper motion can then serve as an anchor point for the interpretation of 
future high-resolution images, which could in principle directly detect the lens or demonstrate 
unequivocally that it is not luminous and so is a sub-stellar object in the foreground disk. 
However, as we now show, such unambiguous results actually require that new images, with 
FWHM~ 50 mas, be obtained "immediately", i.e., before the lens and source have separated 
significantly. 

Suppose, by contrast, the first epoch consisted solely of the 132 mas FWHM images 
already in hand, and that second epoch AO images were obtained 10 years later, which 
(unlike the first epoch) did reach the diffraction limit of 50 mas. The lens will then be 50 
mas from the source, and so separately resolved if it is luminous. But if such a star were 
observed at 50 mas, how could we be certain it was the lens? If the excess light were due 
to an ambient star or to a companion to the lens or the source, then this object could also 
happen to be 50 mas from the source at this latter epoch. For the ambient-star and source- 
companion cases this is obvious. For the lens-companion case, the shear limits discussed 
in Section 18.31 do place some constraints on future lens-companion positions, but as we will 
make clear further below, these still allow it to be 50 mas from the source after 10 years. 

On the other hand, suppose that nothing was detected in such 10-year post-event images. 
In this case, we would know that the lens was not in the bulge, but we would still not be 
able to determine whether the excess light had been due to an ambient star, or companions 
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to the lens or source. And it would be of substantial interest to do so because, in this case, 
a lens companion would be the only clue to the distance (and so mass) of the lens. 

Let us consider now how the situation would change if new images were obtained im- 
mediately with 50 mas resolution, either from HST or using AO. Such images would either 
resolve out the excess flux, or restrict it to 50 mas radius (smaller in the case of HST as 
discussed below). If it were resolved, then the appearance of a "new" star 10 years later 
would have to be due to the lens or its companion. (In principle, such a "new" star could 
be an ambient star that had been hidden at the time of the first epoch, but the probability 
of this is reduced by (50/132)^ = 0.14 and is further reduced by the chance that it would 
happen to be very close to 50 mas from the source at the second epoch.) The proper motion 
of the excess light relative to the source would tell us whether it was an ambient star, a 
companion to the source or lens, which in the last case would give the direction of proper 
motion, thereby confirming that the "new" star was either the lens or a second (and very 
close) companion to the lens. These last two possibilities could not be strictly differentiated. 
However, as discussed in Section 18.31 companions cannot be too close because of the limits 
on shear, so the second companion could potentially be strictly ruled out depending on the 
analysis of the other stars in the image. 

On the other hand, if the first epoch image did not resolve out the excess flux, then, 
as argued above, the ambient star hypothesis would be so much less likely that it could be 
ignored. Appearance of a "new" star in the second epoch would then be either the lens or a 
very close companion to the lens. Again, the strong constraints on the shear would translate 
into very strong constraints on the lens-companion scenario. 

If the first epoch were carried out with HST then these constraints could be tightened 
further. Color- dependent centroid shifts (between say V and /) can be detected for star 
separations down to about 15 mas (assuming an /-band flux ratio of 11%), which (as outlined 
in Section 18. 3p is quite close to the minimum source-lens separation, unless the source is in 
bulge. 

In brief, immediate observations with HST, with foUowup 3 {HST) to 10 (AO) years 
hence, could unambiguously distinguish between a bulge and disk lens, and if the former, 
give a good measurement of the lens mass and distance. Immediate AO observations, if they 
achieved 50 mas, would significantly constrain the possible options, but would not yield an 
absolutely air-tight case. 

We note that the calibrated source magnitude is (V, I, i/)source = (20.76, 19.28, 17.73) ± 
0.05 and the i/-band magnitude of the blend is i/biend = 19-56 ± 0.41. If the blend is in the 
bulge, then (l^,/)biend ~ (24.0,21.9). 
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9.4. Other Bulge Planet Candidates 

The procedures just outlined are challengin g but alternati ve routes to secure detection 



of bulge planets are, if anything, more diffic ult. I G audi I (I20001) discussed the prospects for 
detecting transiting planets in the bulge and ISahu et al. J2006h reported the detection of 16 



candidate bulge planets from a transit survey carried out with the Hubble Space Telescope. 
Two of these were bright enough for radial-velocity follow-up, one of which showed variations 
consistent with a planet with mass m = 10 Mjupiter and the other showed upper limits 
m < 4Mjupitcr- The stars are so bright, however, that their inferred masses indicate that 
at least one (and possibly both) probably lie in the foreground disk. Nevertheless, this 
technique could in principle be pushed harder, particularly when larger telescopes come on 
line. Even then, however, lower-mass planets, m < Mjupiter will probably only be accessible 
with microlensing. 

There are three other planets detected by microlensing for whic h the distances ar e 
neither measured nor st rongly constrained, OGLE-2005-BLG-169Lb JGouM et al.l booeh. 



OGLE-2005-BLG-390Lb JBeaulieu et al.ll2006h . and MOA-2007-BLG-400Lb ( 



Dong et al 



2009tf ). 



In addition, the distance to MOA-2003-BLG-53/OGLE-2003-BLG-235is not precisely enough 
measured to determine unambiguously whether it is in the inner disk or the outer bulge. In 
all four cases, both 6-^ and /i are measured, so we estimate the minimum lens mass that 
would allow the lens-planet system to be in the bulge and the time that must elapse before 
definitive imaging observations can be undertaken. As mentioned below, all of these events 
have large proper motions, fi > 7masyr~^, which generally favor disk lenses. 

For OGLE-2005-BLG-169Lb, 9e = 1.00 ± 0.22 mas and /i = 7 - lOmasyr^^. Even 
adopting the 1 a lower limit on 6'e, then tTj-ci > 75 yuas for stellar hosts with M > Mq. Thus, 
for bulge sources at Ds = 8kpc, the lens distance is no more than 5 kpc. Hence, the lens 
is almost certainly in the disk. Measurements to confirm this relatively secure conclusion 
could be made as early as 7 years after the event, i.e., 2012. 

For MOA-2007-BLG-400Lb, Oe = 0.32 ±0.02 mas and /i = 8.2 ±0.5 mas yr'^ Adopting 
Ds = 8 kpc, the lens would only lie within 2 kpc of the source provided that M > 0.30 Mq. 
Thus, this is a reasonable, but not particularly strong candidate for a bulge lens. The source 
is a moderately bright subgiant, so for 10m class telescopes it is perhaps best to wait for the 
separation to reach 70 mas, which will require about 9 years, i.e., in 2016. 

For OGLE-2005-BLG-390Lb, 9e = 0.21 ± 0.03 mas and /i = 6.8 ± 1.0 mas yr'^. It is 
therefore the best previous candidate for a bulge lens since which is the product of the 
mass and relative parallax, is only a factor 1.6 times larger than for MOA-2008-BLG-310Lb. 
This means that if it were at the bottom of the main- sequence, it would lie about 3 kpc 
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in front of the source and therefore most hkely he in the disk, but if it had significantly 
larger mass it would be in the bulge. However, in this case, the source is a G4 III giant with 
/o = 14.25, which implies Mh ~ —0.85. A lens close to the bottom of the main sequence has 
Mh ~ 11 and so (even accounting for its closer distance) would appear 25000 times fainter 
than the source. While this is an extreme case, it would appear prudent to wait for the lens 
to move 3 FWHM away from the source, which for 10m class telescopes would require about 
20 years, i.e. 2025. If larger telescopes with AO come on line before that, it will of course 
be possible to make the measurement sooner. 

Finally, MOA-2003-BLG-53/OGLE-2003-BLG-235 has a proper motion of 3.3±0.4 mas yr" 
which was sufficient to measure the color-dependent centroid s hift from HST obse rvations 
taken just 1.78 years after the event, but only at the ~ 3 a level ( Bennett et al.|[2006l ). Based 
on this measurement, the lens distance is estimated to be = 5.8l;Q7kpc, which reflects a 
roughly 30% error in the lens-source relative parallax vTrci. Thus this planetary system could 
be in the inner disk or the outer bulge. The color- dependent centroid shift could certainly be 
measured more accurately today, but this would not dramatically decrease the uncertainty 
in Dl, which is fundamentally limited by the 25% uncertainty in 6'| = hMti^qI-, and so in 
TTrei. Thus, pending spectra of this J ~ 21 lens after it is fully separated from the source, it 
will be difficult to prove whether or not this planet is in the bulge. 
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A. Reduction of VLT NACO Images 



Since the reduction of the NACO images is a delicate procedure, we present it in more 
detail. The master darks are median stacked from 5 raw dark frames taken on the same 
night with the same integration time (40 s for H band, 50 s for J and Ks) as the science 
frames. The master fiatfield is obtained from 6 lampfiats taken the same night. A badpixel 
map for correction of the raw frarn es is obtained using the deadpix routine from the ESO 
ECLIPSE package (jPevillardl 119971 ) . The science frames (24 s in J, H and 49 s in Ks) are 
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then dark subtracted, flatfielded, median co-add ed and sky-subtracted using the JITTER 
infrared data reduction software (jDevillardlll999l ). To avoid border effects, we keep only the 
intersection of frames for all the dithered positions for our photometric analysis. 



We use the Starfinder (jPiolaiti et al.ll2000t ) tool to extract the photometry of the reduced 
NACO frames. Starfinder has been specially designed to perform photometry of AO images 
of crowded fields. It creates a numerical PSF template from chosen stars within frame, which 
is then used for PSF-fitting of all stars in the field. Even though the AO correction for the 
given data set is good (strehl ratios of around 10%) and the variation of the PSF-shape 
across the field of view is small, we decide to take star 3 (see Fig. [5]) as PSF template for 
best photometric accuracy on the target, as it is the closest high signal-to- noise ratio star to 
the microlens. 



B. Photometric calibration of IRSF, VLT NACO and H CTIO 

As discussed in Section [TJ there are only two common stars, both with relatively large 
photometric errors, with which to perform a direct photometric alignment between the CTIO 
and NACO systems. As the IRSF images share more common stars with both NACO and 
CTIO, we obtain a more accurate alignment using the indirect transformation NACO-to- 
IRSF-to-CTIO. Specifically, we perform the following steps. First, the IRSF images are 
calibrated with respect to 2MASS reference stars using GAIA/Skycat Fit to obtain initial star 
positions relative to the 2MASS astrometric catalog, and then Tweak is used to refine them. 
We cross identify 1521 objects between the 2MASS and IRSF frames, 779 of which have 
high quality flags (labeled AAA in 2MASS catalog), and then apply two further restrictions: 
keeping only the bright end of the sample, and removing 1.5 sigma outliers . We adopt the 



color terms as given by the IRSF manual and detailed in iKato et al.l (120071 ). and we fit the 
zero point : 

JlRSF.inst = 23.073 ± 0.001 + J2MASS " 0.043( J2MASS - i^2MASs) + 0.018 

i^iRSF,inst = 23.128 ± 0.001 + H2MASS + 0.015( J2MASS - i^2MAss) + 0.024 

i^IRSF,in.t = 22.334 ± 0.001 + i^sMASS + 0.010( J2MASS - i^2MASs) + 0.014 



We apply these relations to the 3006 objects with good cross ID in IRSF images. Up 
to this point, we have calibrated JHK measurements taken by the IRSF telescope. The 
WCS positions of IRSF objects are deduced using the WCSTools routine xy2sky and used 
as references to calibrate NACO images with the WCSTools routine imwcs. In the NACO 
field, we identify 6 bright stars likely not to be affected by blending when comparing IRSF 
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and NACO. Two of them are variable, which leaves us with four stars with the color range 
(J — H) = 0.4 — 0.78. We note that there is no color term in the transformation, and we 
estimate photometric offset between -f^iRSF.caUb and instrumental NACO to be 27.873 ±0.014 
in H. 

We cross-identify 209 stars in the IRSF and CTIO if-band images with matches better 
than 0.8". Wc clip at ±0.1 mag around the mean of Hqtio ~ -f^msF.caiib) and keep 175 
stars. Wc estimate the zero point offset between instrumental i^cTio and i^iRSF.caiib to be 
3.8164 ± 0.0034. 
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Fig. 1.— Top: Light curve of MOA-2008-BLG-310 showing data from MOA (green), /xFUN 
Auckland (orange), /iFUN Bronberg (black), /iFUN SMARTS /-band (red), MiNDSTEp La 
Silla (cyan), and PLANET Canopus (magenta). Also shown is the best fit single-lens model. 
The light curve does not look anomalous at first glance. Middle: Residuals to the best-fit 
single-lens model. Anomalies are apparent at HJD' = 4656.34 and HJD' = 4656.48. The 
noticeable offset in the alignment of Bronberg and CTIO data is an effect of independently 
fitting fs and fb for each observatory (see Section [3]). See Figure [3] for didactic residuals. 
Bottom: Residuals to the best-fit planetary model (the wide solution is chosen for this plot, 
however the close solution is essentially indistinguishable). 
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Fig. 2. — Top: The wide solution Ax^ = 1,4,9 contours in the [ d, q) p lane. The contours 
generated by fixing the limb darkening parameters at the IClaretl (120001 ) values {solid lines) 
are similar to those from allowing F and A to vary freely {dotted lines). In particular, large 
regions of the Ax^ = 1 minima overlap. Bottom: Ax^ contours for the close solution. 
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Fig. 3. — Top: Didactic residuals to the single-lens model. Data points are shown for /iFUN 
Bronberg (black), /iPUN SMARTS (red), and MiNDSTEp La Silla (cyan). The solid lines are 
the best-fit wide (green) and close (blue) planetary models. Bottom: The source trajectory 
(solid black line) showing the extended source (circle) crossing the caustic created by the 
planet at key points in time. The circle radius on the plot is the source radius crossing time, 
t* = pt-E ~ 0.055 days. The caustics for the wide and close models are plotted in green and 
blue, respectively. The density of the caustic points is proportional to the strength of the 
caustic, so that the "solid lines" correspond to stronger magnification while the "dotted lines" 
indicate that the caustic is weaker. The two caustic structures are nearly indistinguishable in 
regions probed by the source. Several of the anomalous features apparent in the the residual 
plot correspond to the limb of the source crossing the caustic. These features are numbered, 
and dashed black lines connect them to the corresponding position of the source. Didactic 
residuals show the difference between the data and a point-lens model that has the same (to, 
Mo, tE, P, /s, fb) as the best-fit planetary model. 
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Fig. 4. — Calibrated color-magnitude diagram of the field containing MOA-2008-BLG-310. 
The clump centroid {circle) is located at [(V^ — /),/]ciump = (1-84, 15.62). The source color 
and magnitude (square) is derived from the best-fit planetary model, [(V^ — /),/]source = 
(1.48 ± 0.01, 19.28 ± 0.05). Assuming the source lies at 0.05 mag behind the Galactic center, 
{V — I)o = 0.69 an d Mj = 3.46, consistent with a post-turnoff G type star, as confirmed 
spectroscopically by lCohen et al.l (120091 ). 
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Fig. 5. — Left: Image taken in H band by the IRSF telescope in South Africa on 2008 
August 4. Right: Median i^-band AO image taken by NACO on VLT on 2008 July 28, 
near the baseline of the event. The lensed source plus blend is indicated as the target and 
reference stars are circled and numbered. 
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Fig. 6. — Probability distributions (normalized to unity) for the model-derived magnified 
source flux at the time of the NACO image. The curve for the case of no parallax is plotted 
in black and unconstrained parallax is plotted in gray. The mean for each distribution is 
indicated by a dotted line, and the standard deviation in each case is 0.04 mag. The best 
estimate of the target flux on the NACO image (21.28 mag) is marked by the vertical black 
line, and dashed lines are 0.05 mag conservative error bars. The error bar at 21.33 mag can 
be considered a robust lower limit on the amount of light detected on the NACO image. 
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Fig. 7. — Search for additional shear due to distant third body, parametrized by the angle 
between the planet and the third body, 0, and the width of the resulting Chang-Refsdal (1979, 
1984) caustic: Wcom = 4gcom'^ram- Open symbols represent improvement relative to xLin ~ 2, 
i.e. the value expected from adding two degrees of freedom to the single-planet model. Filled 
symbols indicate worse x^- Colors (black, red, yellow, green, cyan, blue, magenta,white) = 
(1, 4, 9, 16, 25, 49, > 49). There is a weak detection of shear at (logiycom, 0) = (—3.28,40°), and 
a 3cr upper limit logWcom < —3 over most angles. Note that the angle between the shear 
direction and the direction of source motion is + a, where a ~ 69° is the angle between 
the planet direction and the direction of source motion. 
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Fig. 8. — Mass versus equilibrium temperature for planets detected via radial velocity (cir- 
cles), transits (triangles), imaging (stars), astrometry (squares), and microlensing (hexagons). 
If the blended light aligned with the event that was identified by VLT NACO is in fact due 
to the lens, MOA-2008-BLG310Lb would be the first microlensing detection to fall on the 
Snow Line. The tail on the marker for this detection indicates where the planet might fall 
if the host is a lower-mass star, rather than being identified with the blended light. (Data 
taken from ,http: //exoplanet.eu/^ maintained by J. Schneider.) 



